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INFRARED DIFFERENTIAL ABSORPTION FOR ATMOSPHERIC POIiUTANT DETECTION 
I. INTRODUCTION 

This first semi-annual report summarizes the progress made to date 
in the generation of tunable infrared radiation and its application to 
remote pollutant detection by the differential absorption method. 

At the Laser Radar Conference in Sendai, Japan held recently, it was 
recognized that future remote pollutant measurements depended critically 
on the availability of high energy tunable transmitters. Furthermore, due 
to eye safety requirements, the transmitted frequency must lie in the 1.4 ^im 
to 13 ^m infrared spectral range. 

Recent experiments by Walther et al. , at the University of Cologne in 
1 2 

Germany and by Grant et al. , at S.R.I. have shown experimentally that the 
differential absorption method offers the best potential for sensitive remote 
pollutant detection. There was general agreement at the Sendai conference that 
this is the case. 

Our paper, describing progress in our high energy tunable coherent source 

was the only paper presented at the conference concerned with the generation of 

high energy tunable infrared radiation. A preprint of the paper describing 

the 1.4 \m to 4.4 LiKbO^ parametric oscillator is given in Appendix I. 

4 

Appendix II is a preprint of a review paper on remote pollutant detection 
written for the ICO Meeting held in Tokyo just prior to the Sendai Meeting. 

This paper has been submitted to Opto-Electronics for publication. 
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II. DIFFERENTIAL ABSORPTIOM METHOD 


We have given further consideration to pollutant measurement by the 

differential absorption method. Three areas of concern have been identified. 

They are required transmitted energy vs wavelength in the infrared due to the 
2 

approximately X decrease in Mie backscatter with wavelengths between 1 jini 
and 10 pm; the decrease in scintillation in the infrared and the reduction in 
signal to noise ratio due to digitatioh noise. 

2 

The decrease in returned signal due to the X decrease in Mie scattering 

-is being investigated in more detail theoretically. But at this time differential 
absorption measurements may be limited to wavelengths less than 5 with the 
present transmitter operating at 20 mj pulse energy. 

A review of scintillation studies seems to indicate the likely possibility 
of an alternate on and off resonance pulsed system as opposed to a double pulsed 
system. In the infrared^ scintillation decreases in amplitude with increasing 
wavelength. Due to the high cost and difficult technology involved in a double 
pulse system we plan to carry out initial work with the alternate pulsed approach. 

Recent advances in A to D converters has led to a 10 bit, 10 MHz system 
at a cost of under $10,000.^ This A to D converter, followed by a buffer 
memory Is fully capable of processing the return radar signal in conjunction 
with our PDPllE-10 computer. We feel that the 10 bit capability is essential 
if full sensitivity is to be obtained in the differential absorption measurements. 

We are presently discussing possible additional support with the Electric 
Power Research Institute (EPRi) to aid in the purchase of the A to D converter 
for the data receiving system. 
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HI. HIGH ENERGY WIDELY TUNABLE INFRARED SOURCE 

We have successfully demonstrated the 1 . 14 - p.m to h.k LiNbO^ parametric 
oscillator source. The important results are summarized in the preprint 
presented in Appendix I. 

Our present effort lies in the evaluation of the 1 . 0641 m amplifier 
chain and operation of the oscillator at higher input energies. We expect 
to complete this work within the next two weeks. 

Our schedule calls for the completion of the computer controlled oscillator 
stage by January 1975* We should be ready to initiate external pollutant 
measurements by early Spring. 

The growth of (01*^4) LiNbO^ crystals has been carried out by three 
companies so that a commercial supply of these crystals should be available soon. 

XV, DATA PROCESSING SYSTEM 

The PDPllE-iO minicomputer was received In July. It is now fully operational 
with KT-11 and basic software. Recently we received the CAMAC interface crate 
and stepper motor controllers. A software program was completed which allows 
programing of stepper motors and D to A converters. We have demonstrated 
the capability by programing a stepper motor. 

We are presently evaluating data display systems. Costs ‘and capabilities 
range from $3,500 for a CAMAC controlled color TV display to over $25,000 for 
a full computer controlled random access color display. 

At the Sendai conference C. Frush of NOAA presented a paper describing a 
computer compatable fast A to D converter with gain switching end buffer 
memory. We have called Dr. Frush and ore more carefully evaluating his processor. 
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At this time, the specifications of Computer Labs. 10 bit A to D converter 
appear to be the best available. Vfe hope to complete the evaluation of the 
available systems during the next two months. 

V- MEASUPEMENT PROGRAM 

We have initiated design of a telescope receiving system and the 
necessary building modifications to install a receiver-transmitter system 
on the Microwave Laboratory roof. From our roof location probable targets for 
a future air pollution study are down-town Palo Alto (range ~lkm), Stanford 
parking lot (range ^.-200 m) and the Stanford Steam Plant (range m). We 

also have numerous targets for topographical reflectors including Hoover tower 
(range 600 m). Eucalyptus trees (range 500 m) and the hills across the Bay 
(range *^10 km). We will more fully evaluate these potential targets as our 
measurement program begins. However, at this time the natural gas burning 
Stanford Steam Plant appears to be an ideal active source for future study. 
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APPENDIX I 


A 1.4 ^l^l to 4.4 \m HIGH ENERGY ANGLE TUNED LiNbO 
PARAMETRIC OSCILLATOR 

E.L. Herbst, R.N. Fleming and R.L. Byer 


ABSTRACT 

We have operated a high gain, angle timed> singly resonant LiNbO^ 
parametric oscillator pumped directly at 1. 06 pm by a switched Nd:YAG- 
laser. The oscillator angle tunes from degeneracy and operates over the 
entire 1.4 pm to 4.4 pm range., Output energies of greater than 1 mj/pulse 
at 5 pps have been observed with a 15jg energy conversion efficiency. The 
key-to this device is the large LiNbO^ crystals fabricated from new [01-4] 
grown boules. Crystals up to 15 mm diameter and 5 cm in length have been 
cut at the nominal 4?° orientation. These crystals will potentially handle 
over 2J of optical energy at 1*06 pm. 



A 1,4 to 4.4 tun HIGH ENEEGY I^GtK TUNED LlNbO. 
PAEAMETRIC OSCILLATOR 


At present two methods are available for obtaining "peak powers greater 
than several hundred watts of continuously tunable coherent radiation in the 
X.5 to 4.5 jutt region. These are parametric mixing and parametric oscillation. 

Since in parametric mitring two input wavelengths are needed one of which is 
tunable^ mixing wavelengths have been limited primarily to the visible and 
near infrared where dye lasers provide the tunable component. Although high 
peak input powers are available efficient down conversion has not been realized 
because of problems of beam overlap and phase mismatch. * * parametric 
oscillators, potentially efficient converters in the region, have been 
limited to low power or low efficiencies because of operation far from 
degeneracy, small crystal size or low ciystal damage threshold. 

We have recently studied a potentially high energy widely tunable infrared 
source* The device is an angle txined, singly resonant LiNbO^ parametric 
oscillator pumped by a Q-switched I .06 jim Kd:YAG laser. The oscillator operates 
at room temperature and txines the complete range of 1.4 nm to 4.45 pja. The 
key element, in. this device is the large LiUbO^ crystals available by using a 
new [01*4] boule growth direction. r“ 

LiHbO^ boules grown along the [01*43 axis^^ yield large diameter, strain 
free, high optical quality crystals. The boule axis lies 38 ° to the optic axis 
in the yz plane. Thus large parametric oscillator crystals can be fabricated 
at the nominal 47° phasematching angle desired for a I .06 pm pumped parametric 
oscillator. To date we have grown fourteen boules and have fabricated oscillator 
crystals I .5 cm diameter by 5 cm long. Using SIO^ anti-reflection coatings, these 
liNbO^ crystals have operated in the parametric oscillator at greater than 

o 

2 j/cm energy density without damage. For our 20 nsec Q- switched pulse the 


7 


, g. 

corresponding surface intensity is 100 MW/cm • 

The boules are grown by the Czochralski method at the congruent melt 

exposition using seeds with a known [01*0] direction. After annealing, 

the boules are poled to maximize the effective nonlinear coefficient. Since 

^eff * ^1 + P) + dgg cos(f3 + p) sin 3cp where ® is the phasematching 

angle, p the double refraction angle and cp the orientation with respect 

to the [11-0] direction, and d^^ and are known to be of opposite 

aign, the effective nonlinear coefficient is maximized for propagation 

•12 

in the second and fourth yz quadrants. For [ = 6.25 x 10 (m/v) 


.-12 


and Idggl = 3-3 X 10 *“ m/v , = 6-88 X lo"^ m/v for 6 = 4?” and 

p = 0.04 rad- Thus the effective nonlinear coefficient is slightly larger 
at 6 = than that for 6 = 90° . 

The experimental setup, is shown in Fig. la. The pump source is a low 
energy flashlamp pumped Nd:TAG laser operating at 1.06^ at 10 pps. The laser 
Is ^ectro-optic Q- switched with output energies of up to 20 aiJ In a 20 ns pulse. 
Operation is limited to a TEM^mode with an internal aperture but no effort was 
made to limit the spectral bandwidth. Figure lb shows a photograph of the , 
-llHbO^ oscillator with 1 cm. diameter by 5 cm crystal within the oscillator cavity 
The laser output having a 76O ti spot size is incident directly on the 
oscillator cavity without collimation or mode matching. For operation over 
the entire tuning range two sets of parametric oscillator mirrors were used. 
However, broadband mirror coatings covering the 1.4- to 2.1 ^.m region will allow 
tuning over the entire range from 1.4 ^im to 4-. 4- jam with one set of optics. 

The oscillator cavity consists of two long radius or flat mirrors spaced by 
7 cm. The LiNbO^ oscillator crystal was held in a copper block at room 
temperature. No attempt was made to stabilize the crystal temperature and 
variations in oscillator wavelength with changes in room temperature 


0 



occurred. For these initial experiments these variations were not a problem 

o 

hut for long term stability an oven stabilized to j; .1 C will be required. 

Figure 2a shows the calculated timing curve for the LillbO^ parametric 

oscillator using Hobdens and Warner ^s index of refraction relations. Also 

plotted are the measured data points. The triangular points were taken with 

a 2 m high reflector input mirror and a flat output coupler. These were 

standard mirror coatings centered at 1.8 ^im. The single circle data point 

at \m was taken with standard 1.32 \im Wd:YAG laser optics consisting of 

a high reflector flat input mirror and a flat 10^ output coupler. The 

oscillator tuned continuously over the entire range. Tuning of the oscillator 

beyond the pm point was limited by the inability to reach threshold due 

to Increasing infrared absorption of the idler wave. The oscillator threshold, 

measured near 1-9 pm. was 4 mj corresponding to 1-5 X ^0^ watts of pump power. 

-1 • -1 

The measured gain bandwidth varied from 3 cm near I .3 to 23 cm near 
degeneracy. 

Because of' the large double refraction in LiWbO^ at the li-7® crystal angle, 
8 large pump beam spot size is needed to reduce the effect of pump beam walk- 
off on parametric gain. An estimate of the pump beam spot size can be obtained 
using 

1 



( 1 ) 


where p is the double refraction angle. For a 5 cm LiWbO^ crystal at 47^ 
with p= 0.04 radians, w^ = 8 OO p . Using higher power lasers and maintaining 
constant intensity the spot size increases such that the parametric gain 
approaches the plane wave limit given by^^ 



vfaere 




W3^ 


L -Zi 

=3 A. 


( 3 ) 


and i is the crystal length and = it w^/2 . At degeneracy we have 


s 2.3 X 10“^ 

^3 


Q 


cm 


-2 


Thus for F,pr-’"g near the damage intensity threshold of 80 W/cm with large 
spot size -beams, single i®ss gains of G = 3./^^ e^^^ = X 1°^ are possible 
In 5 cm ca^stals. These gains are approximately that required for superradiant 
operation and illustrate the high gains possible in par^etric 'ievices. 

An important aspect of using a large spot size to pump a parametric 
oscillator that is often overlooked is the ability to use tilted etalons 
inside tlie oscillator cavity with low insertion loss. For a tilted etalon 


the insertion loss is given by 


16 


R 


Aod^ 


(1 - R)^ \nw / 


• (h) 


where a is the tilt angle, d the etalon thickness, R the etalon reflectivity, 
and w the beam radius. The requirement for large spot sizes at the resonated 
wave is apparent. Figure 2b shows the spectral output of the oscillator with 
a 2 mm thick, finesse of ten etalon inside the oscillator cavity. The line- 
width of the individual etalon modes is less than 0.1 cm ^ and is not resolved 
by the spectrometer. By scanning the etalon angle we tuned an individual mode 
over many free spectral ranges without increasing the cavity loss. The amplitude 
envelope of the etalon modes is determined by gain lipewidth of the parametric 

oscillator. 
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Figure 3 demonstrates the oscillator output stability. Figure 3a * 
shovs that the peah-to-peak fluctuation of the total oscillator output 
over a five minute period is less than 3«g as monitored by an Eppley 
thermopile with a one second time constant. During this measurement the 
oscillator operated at 5 Pps with 1 mj per pulse output energy. Figure yo 
shows 1000 consecutive output pulses at 20 pps. The oscillator operated 
7 fflW average power with a peak to peak stability of 3 ^. Since the InSb 
detector did not resolve the 20 nsec pulse, the oscilloscope trace ‘is a 
measure of output energy stability with a 1 jisec time constant. 

In. conclusion, we have demonstrated a widely tunable infrared source 
that can be rapidly and conveniently angle tuned over the entire 1.4 pm to 
4.4 |jua range. It is significant that the source requires only one set of 
optics and a single LiBbO^ crystal to tune the entire range. Conversion 
efficiencies of 1^ and operating stability of better than 3 ^ have been 
observed. The availability of large I 5 mm diameter by 5 cm long, 47° cut 
LiKbO^ crystals allows operation of the oscillator at up to 2 J of I .06 pm input 
pulse energy. The large oscillator mode radius allows the use of a tilted 
' etalon for line narrowing without significant reduction of output power. 

Finally, it should be noted that the oscillator covers more than one 
octave of tuning range so that its frequency can be extended. into the visible 
and ultraviolet by second harmonic generation in LillbO^ and LilO^ and still 
maintain complete spectral coverage. Frequency extension into the infrared 
is most easily done by mixing the signal and idler waves in semiconductor non- 

17 18 

linear crystals such as AgGaSe^ * and CdiSe which phasematch over a 3 [am to 
IS pm and 10 pm to 25 pm range. The high output power of the LiBbO^ tunable 
source insures that the nonlinear doubling and. mixing steps will proceed with 
good efficiency. 
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FIGUEE CAPTIONS ‘ 

« 

X* a) Schematic of experimental arrangement showing the electro- optic 
Q-Swltched 1.06 \m Nd: YAG laser and LiNbO^ parametric oscillator, 
b) Jhotograph of the singly resonant LiNbO^ parametric oscillator 
vitb a 1 cm by 5 cm crystal within the oscillator cavity. 

2* a) Calculated parametric oscillator tuning curve and accurately 
measured wavelength .vs crystal angle points. The entire tuning 
range to the circles was continuously scanned using a single 
15 BEQ X 5 cm oscillator crystal. 

b) Spectral output of the oscillator with a tilted etalon. The 
llnewidth of individual etalon modes is less than 0.1 cm ^ and is 
not resolved by the spectrometer. The amplitude of the etalon modes 
shows the gain envelope of the parametric oscillator. 

‘3; a)- Average output stability of the oscillator operating at 5 pps, 

1 mJ/pulse for five minutes with an Eppley thermopile, 

b) Peak to peak stability showing 1000 pulses at 20 pps. The average 

and peak output fluctuations are less than 
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REMCTTE POLLUTION MEASUEEI-IEKT 
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♦ * * 

ABSTRACT 

The availability of high energy tunable laser sources has extended 
the possibilities for remote air pollution measurement. This paper presents 
a discussion and comparison of the Raman method, the resonance and fluorescence 
hachscatter method, long path absorption methods and the differential absorption 
method. A comparison of the above remote detection methods shows that the 
absorption methods offer the most sensitivity at the least required transmitted 
energy. Topographical absorption provides the advantage of a single ended 
measurement, and differential absorption offers the additional advantage of a 
fuHy depth resolved absorption measurement. Recent experimental results 
confirming the range and sensitivity of the methods are presented. 

Future widespread use of remote pollution monitoring will probably be 
limited to the infrared to interact with molecxilar spectral bands and to 
meet eye safety requirements. "" 



REMOTE AIR POLLUTION MEASyEE2»EKT 

!• IHTRODUCTION 

A number of review articles on air pollution measurements have already 
been published. However, recent advances in tunable laser technology 
coupled with experimental verification of remote sensing methods bring remote 
pollution measurements very close to a practical state. This paper reviews 
the field of remote pollution measurement and indicates the requirements for . 
tunable laser transmitters. 

The adverse effects of polluted air on human health is well established 

and an estimation of the cost savings for human health alone in the United 

T 

States is over one billion dollars for a 5^ pollutant reduction. The 

KnVironmental Protection Agency has established maximum pollutant levels and 

goals for future reduction In pollutant levels. To meet these standards, 

accurate sensitive measurement methods must be established. Major pollutants 

of interest Include ozone, nitrogen oxides, carbon monoxide, sulfer dioxide 

and hydrocarbon compounds. These pollutant molecules are present in urban 

8 

atmospheres in concentrations from .01 ppm to over 10 ppm. Remote pollutant 
measurement methods must have measurement sensitivities in this range If they . 

are to be useful for atmospheric monitoring. In addition, the detection methods 

« 

must be capable of monitoring molecular vapors. Both Raman scattering and 
infrared absorption on molecular vibration bands satisfy this last requirement. 
The detection sensitivity inherent in absorption and the regard for eye safety 
requirements probably limit future remote pollution detection methods to the 
infrared spectral region. 
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A. Historical Overview 

Historically lasers were first used in radar type remote sensing 

applications. The acronyn LIDAR for Light Detection and Ranging was 

9 

applied to laser ranging systems. In 19^3 Fiocco and Smullin recorded 

backscattered echoes from turbidity in the upper atmosphere and Ligda^^ 

measured- backscatter from the molecular atmosphere and haze. Since that 

early work considerable effort has gone into LIDAR development and applications. 

The remote pollutant detection methods considered in this paper are an extension 

of LIDAR and therefore work in LIDAR and its applications are of interest. 

Early reviews and discussions of LIDAR are given by Inaba^ Kobayashi and 

lehimura ^ Collis , Kent and Wright and Collis. More recent reviews 

including discussions of the wide ranging applications of LIDAR are given by 

DeCT and Little^. , Strauch and Cohen^^ , Collis and Uthe^^ , Allen and Evans^*^ , 
^ ^ 19 

Hall and Collis and Russel. The field has grown rapidly and has laid much 
of the ground work for remote pollution detection methods. 

Using the then available fixed frequency laser sources Raman scattering 

20 

was first used to detect molecular constituents of the atmosphere by Leonard ^ 

21 22 2 ? 

Cooney and Melfi et al.. Later Inaba and Kobayshi ^ discussed pollutant 

detection by Raman scattering and presented experimental work showing the 

24 

detection of SO^ and CO^ molecules. Remote detection by Raman scattering 
has' the advantage of using a fixed frequency laser source which can operate 
in the visible or violet spectral region where photomultiplier detectors are 
available. Unfortunately^ the molecular detection sensitivity possible by 
remote Raman scattering is not adequate for trace contaminant measurements. 
However, considerable effort has 'gone into using Rrman scattering for monitoring 
water vapor profiles^^ and major components of the atmosphere. Recently 


- 23 - 


T 



Inaba and Kobayshi^ reviewed Laser-Raman Rrdar and its applications to reniote 
monitoring. Remote detection by Raman backscattering is discussed in Section 
II. 

The availability of tunable laser sources opened additional remote pollutant 
detection possibilities.^ They include long path absorption^ detection of 
resonance fluorescence and differential absorption. The advantages of these 
detection methods is the increased molecular detection sensitivity due to the 
large absorption and fluorescence cross sections. In addition^ absorption and 
fluorescence methods can be used for the detection of both molecular and atomic 
species. 

28 

Using multiple line laser sources Hanst and Morreal and Snowman and 
GiUmeister^^ investigated the possibility of pollutant detection by long path 
ab^rption. The low required transmitted powers for long path absorption also 
allows the use of tunable cw infrared diode laser sources-^ and Raman spin flip 
laser sources, In addition detection sensitivity can be improved by heterodyne 
techniques. ^^Although long path absorption provides the most sensitivity of any 
remote detection method it has the disadvantages of being double ended and lacking 
depth resolution. The first disadvantage can be overcome by obtaining a return 
signal from topographical targets but at an increase in required transmitter 

power. Both disadvantages are overcome by the differential absorption method 

02 5 

first discussed by Schotland and analyzed in detail by Measures and Pilon 

33 

and Byer and Garbuny. 

The differential absorption method provides a fully depth resolved absorption 
measurement of pollutant concentration vs range by using the molecular and Mie 
particle scattering as a distributed reflector. The required transmitter power. 



£tlgaaX to noise ratio of the returned signal^ and sensitivity of the differential 
abosrptiOQ method is discussed in detail in Section II. 

The improvements in tunable laser scoirces has proceeded rapidly so that 
tunable sources are now considered available for remote pollution measurements. 
Without the use of tunable laser sources, remote air pollution detection is 
limited to LIIiAR and Laser-Raman methods neither of which is adequate for the 
measurement of weahly concentrated molecular pollutants. Tunable laser sources 

increase remote pollutant measurement possibilities to include resonance fluor- 
escence and the absorption methods which do provide the capability of sensitive 
remote detection- The importance of tunable laser sources is evident in 

discussion of remote detection methods presented in Section II. 

Table I lists remote pollution detection methods and provides an outline 

of the measurement methods discussed in Section II. 

B- Atmospheric Parameters 

The atmosphere attenuates a transmitted optical beam by elastic and in- 
elastic scattering and by absorption. Some of the scattered light returns to 
a receiver to provide a signal from which atmospheric properties can be determined. 
In general an Incident beam of intensity I^ propagating through the atmosphere 
is attenuated exponentially or 


1 = ^ 0 ® 




( 1 ) 


vhere I is the transmitted Intensity over a distance r and is the 

atmospheric volume extinction coefficient which is composed of a sum of terms 


'“mIE “aBS 


( 2 ) 
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TABLE I 


REMOTE AIR POLLUTION i^EASUREMENT METHODS 


ELASTIC SCATTERING 

Rayleigh 

me 

INELASTIC SCATTERING 

Raman 

Resonance fluorescence 
Fluore scence 

ABSORPTION 

Long path - retroreflector 
Long path - topographic target 
Differential absorption 
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0 ^^ and 0:^2 are the Rayleigh, Mie and absorption extinction 
coefficients. The Rayleigh extinction coefficient is due to molecular scattering 
and can be calculated as a function of altitude for known atmospheric molecular 
densities, Mie scattering is due to particulates and aerosols in the atmosphere 
and is highly variable in both particulate size, wavelength and particle distrib- 
ution. For identical molecules of density N , the volume extinction coefficient 
can be written in terms of a cross section a = oN where a is the total cross 
section (cm ) for the optical interaction with the molecule or particle. 

In general the scattering of light from molecular or particulate matter is 

not isotropic. The cross section is then specified by the differential cross 

dC 

section a(9, p) =: — where 9 is the scattering angle and cd is the polar- 
ization angle. For isotropic scattering the scattering cross section per unit 
solid angle is da/(^ s o/kit • For backscattering (9 » tt) the volume backscatter 
coefficient 0 = Hh(jt) . 

For molecular Rayleigh scattering the differential cross section is"^ 



^ 2 , 2 ,.2 - 

« (n “ 1) ^ 2 2 2 

— & _ — (cos 9 cos p + sin p) ( 3 ) 

dn 3 X. IT 

where n is the refractive index, N the density and \ the wavelength. For - 

o 

unpolarized light the angular factor averaged over p reduces to (cos 9 + l)/2. 
Finally, for backscattering (9 = jt) the Rayleigh backscatter coefficient is 
uniquely related to the volisne extinction coefficient by 


so that it is possible to determine N by Rayleigh backscattering. 

In the above discussion the Raman volume extinction coefficient has been 
neglected since it is at least an order of magnitude less than the Rayleigh 
extinction coefficient. However, for reference and later use, the Raman 
differential cross section on resonance is given 


^RAM ^ ^ ^RAM 

dD 3t 


(5) 


where and oig are the laser and Stokes frequencies, ^rAM 

half width of the Raman line and is the imaginary part of the Raman 

susceptibility. A list of Raman frequency shifts and cross sections is given 
by Inaba and Kobayshi. The differential Raman scattering cross section is 
approximately 10 ^ cm /sr for diatomic molecules. 

The large variation in size, density distribution and properties of 
atmospheric aerosols preclude an accurate calculation of the Mie scattering 
cross section. However, because of its importance in the differential absorption 
method, it is useful to summarize the properties of light scattering from 
atmospheric particulates and aerosols. Scattering from particulates is 
characterized by the size parameter a 2jtr/\ where r is the particle radius. 
The ratio of scattering cross section a(n,a) to geometric cross section is 




a^(n,a) 


It r 


For a uniform concentration of spherical particles the volume extinction 

2 * 

coefficient is = iir k^(n,a)N> . The factor can be calculated in 
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tferms of an infinite series as first obtained by Mie. Mie 

scattering has been treated by van de Hulst,^^ Kerker^^ and Born and 
1^0 

Wolf^ and computer calculations have been used to investigate the detail 
properties of Mie scattering. An example of the rapidly varying scattering 
amplitude as a function of a is given by Strauch and Cohen. 

In the atmosphere the particle sizes are distributed and a functional 
dependence for the density dependence on radius N(r) is needed to evaluate 
the total scattering. Junge^^ proposed the distribution dlf(r) = C 
for aerosol distributions in the atmosphere with v^ 3 . Elterman^^ has also 
evaluated the distribution of atmospheric aerosols. Fortunately^ the complicated 
integrals over particle size and angular parameters necessary to calculate Mie 
scattering cross sections^ can be avoided by the use of empirical relations. 

For example, the visibility and average extinction coefficient obey the 
rel«ition^^^**^^ 

V const/ p 


vhere V is the visibility and 0 is the average extinction coefficient over 
the range. 


A useful empirical relation relating the Mie backscatter coefficient to 

kg 

visibility is given by 



idiere \ is in microns and V is in kilometers. This expression is an 

approximation and may not hold under unusual atmospheric conditions. However, 

hi U8 

recent measurements ' have shown that Eq. (6) does not apply to Mie scattering 
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with laser sources but with increased variation in backscattered signal 

compared to white light sources. 

Extension of Eq- (6) into the infrared is questionable since the particulate 

size distribution decreases rapidly for particle sizes greater than 1 to 2 iim. 

However, the Mie backscatter coefficient has been recently calculated for 

49 50 

aerosols, fog and rain conditions for wavelengths out to 10 jam. The 

-2 

calculations show that varies approximately as X in the infrared 

-4 -1 

between 1 jim and 10 jim reaching 10 - 10 ^ km at 10.6 |im. Recent measurements 

51 -4 “1 

using a COg laser source'^ show that 10 km at 10.6 |im confirming the 

calculated values. For a more complete discussion of atmospheric scattering see 

52 53 

Goody and Kondratyev. 

The final parameter in the atmospheric volume extinction coefficient is 

absorption. In the ultraviolet spectral region atmospheric absorption is not 

o o 

serious for wavelengths longer than 2500 A. Below 2500 A absorption due to 

atmospheric oxygen becomes important with an absorption coefficient reaching 
“1 o 

1 km at 2450 A. The scattering contribution to the volume extinction 
coefficient in this region is also significant being between 1 km"^ and 5 km'^ 
for visibilities of 10 km to 1 km.^ 

In the infrared absorption is the major contribution to the volume 
extinction coefficient. The atmospheric transmittance in the infrared has ^ 
been subject of a number of investigations. Within the infrared atmospheric 
windows the absorption loss over a 1 km pathlength is significant even under 
clear day conditions. For example, for a 5 km visibility, a relative humidity 
of 755? Slid sn air temperature of 60*^F, the absorption coefficients in the two 
infrared windows 2.7O to 4.20 ^m and 4.40 to 6 [ixa are approximately 0.4 and 
0.8 km The scattering coefficients are almost an order of magnitude less. 
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Trom this discussion we conclude that the atmospheric volume extinction 

coefficient is approximately 1-2 km ^ in the ultraviolet and visible and 

0.5’ tan ^ in the infrared window regions. The magnitude of the extinction 

coefficient sets a maximum remote sensing range for a horizontal path of less 

than 10 km if significant transmitted beam attenuation is to be avoided. 

Two other properties of the atmosphere are important in remote sensing. 

They are atmospheric scintillation and background signal levels during both 

daytime and nightime operating conditions. A large amount of work has gone 

toward the characterization of atmospheric spatial scintillation due to 

57 

turbulence. For the propagation of laser beams through the atmosphere 

spatial scintillation leads to beam distortion, displacement and angular 

variations. Hansen and Madhu"^ have considered laser beams reflected from 

a tetroreflector target and have shown that atmospheric phase fluctuations at 

the retroreflector lead to angle scintillations of the beam. Dabberdt and 

59 

Johnson have recently studied the wavelength dependence and saturation 

characteristics of spatial scintillation.^^ As expected atmospheric scintillation 

is less for infrared than for visible wavelengths.^^ 

Temporal scintillation of the atmosphere is also important in remote 

55 

sensing. Previous work'^-'^ shows that temporal scintillation is constant to 
near 20 Hz and then falls linearly to zero at near 200 Hz. For remote ranging ” 
with nanosecond pulses over microsecond round trip pulse range times the 
atmospher is essentially "frozen" and does not lead to pulse amplitude variations. 
For cw measurements or slow pulse repetition rate experiments, pulses must be 

averaged to reduce the 1-6^ scintillation modulation of the transmitted beam.^^ 

# 

Increasing the receiver aperture size also reduces the scintillation induced 
modulation. 
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Atmospheric background radiation is of concern in evaluation detector 
signal to noise and system receiving parameters. Background radiation includes 
elastic and inelastic backscattered radiation as well as natural background 
radiation. Here we assume that background signals from Rayleigh, Raman, Mie 
and fluorescence scattering can be estimated from previous discussions and 
consider only that due to natural sources. 

The background power at the receiving mirror, due to an extended source 
filling the field of view, is^*^ 

\ \ (7) 

^ere T^ is the atmospheric transmittance at \ the optical bandwidth 
of the detection system, the receiving mirror field of view and B(\) 
the spectral radiance of the' background source. The spectral radiance of the 
clear daytime sky peaks in the visible due to scattered solar light at a value 
Of approximately, lo'^ Wcm’^ . In the infrared the spectral radiance 

again peaks due to thermal radiation near I3 at approximately 10*3 wcm*^ 
pm sr . Below 300O X the osone absorption layer in the upper atmosphere 
screens the solar radiation thus producing an effective night condition. 

Eq. (7) shows that the background signal is reduced by decreasing the field 
of view, the acceptance bandwidth and the receiver area. For a shot noise 
limited detector, such as a photomultiplier, the background power must be kept 
less than the signal power to maintain sensitivity. However, for dark current 
limited detectors such as infrared photodetectors, the background power incident 
on the detector can be greater than the signal power without degrading the 
sensitivity. The upper limit on the received bac.kground power is determined 
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by background fluctuations or detector saturation. The filtering of unwanted 
background radiation is thus much less of a problem for infrared remote sensing 
than for visible sensing. 

The final parameter of concern in remote sensing in the atmosphere is 

the eye safety problem due to the transmission of high peak power laser beams. 

Recently such concern has been quantified by the adoption of eye safety standards 

by the American National Standards Institute. The reccmmended eye safety 

Maximum Permissible Exposure (MPE) in the visible band, O.h ^im to 0.694 jam is 

5 X 10 j/cm j increasing linearly to 2.9 X 10 j/cm at 1.06 um and maintaining 

that level out to 1.4 ^m. In the infrared range 1.4 ^im to I3 |xm the MPE is 

0.56t^/^ j/cm^, which yields 14 X lo'^ j/cm^ for a pulse length t of 300 nsec. 

-4 2 ^ 

In the ultravoilet the MPE varies from 2 X 10 j/cm at O.3IO jam to 3 X lo” 

2 ■*4 2 

j/cm at 0-26 jam to 10 j/cm at 0-20 pm. Between 0,31 pun and 0.40 pirn the 
MPE is 10"^ j/ca?. 

_ The very low MPE in the visible region rules out the use of high energy 
pulse laser sources for remote measurement. Since all molecular absorption 
lines in the infrared lie beyond 1.4 jim, eye safety standards can be easily met 
ly use of a small beam expanding telescope. Public safety and eye safety 
standards will become increasingly important as remote sensing measurements move 
from the research stage to practical applications in the open environment. ^ 

n. MEASURE^^ENT METHODS - A COMPARATIVE ANALYSIS 
A. General Formulation 

The remote pollutant measurement methods are listed in Table ,1. In this 
section a general formulation is presented for the analysis of the sensitivity, 
range, and S/N for the pollutant detection methods. 
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The backscattered signal from a target at range R is given by 



•Hhere P and P- are the received and transcaitted powers, K is the optical 
r o 

system efficiency, A the area of the receiving telescope and (p/:t) the 
effective reflectivity of the remote target. The transmitted radiation inter- 
acts with the pollutant molecules through the atmospheric volume extinction 

coefficient Cf defined in Eq. (2). Returned signals due to Raman, resonant 

A 

fluorescence and fluorescence are included in the effective target reflectivity 
. The strength of the returned signals by the inelastic scattering processes 
is related to the scattering cross section for the process. The depth of the 
atmosphere probed by backscatter methods is shown schematically in Fig. 1. The 
campled depth AR is given by 

*E = (9) 

where r - t + t , + is the sum of the laser pulse width, j > "the 
• p d F P 

detector integration time, molecular fluorescence time r-p - 

For Raman, ultraviolet fluorescence and differential absorption measurements — 
is negligible and the sampled depth is set by either the laser pulse width 
or the detector integration time. Of course the particular sampled depth at 
range R is determined by the standard Radar time of flight measurement. 

Table H lists the scattering and absorption cross sections for the 
measurement processes of interest*. It is readily apparent that absorption 
and fluorescent cross sections are about ten orders of magnitude larger than 



r 



TABLE II 


VJ1 


EEFRESENTATIVE CEOSS SECTIONS 


Process 

Cross Section 

Lifetime 

Scatterer 

Spectral Range 


da . 2. 

(sec) 


(pm) 

ELASTIC SCATTERING 





Rayleigh 

Mie 

10-^7 

10-21 . 10-2‘‘ 


' °2 

particulates, 

aerosols 

0.20 w 1.0 
0.20 « 10 

INELASTIC SCATTERING 





Raman 

10-3° 


all molecules 

uv, visible 

Resonance 
fluorescence and 
Fluorescence 

' 




uv, visible 
IR 

lo"^*^ - lo"^"^ 
10-18 - 10-23 

10-8 - 10-8 
lo"i - io"8 

atoms, molecule s 
molecules 

uv, visible 
2-20 

ABSORPTION 





uv, visible 

IR 

lo"^^ - 10*^*^- 
10‘ lo"^^ 

1 


atoms, molecules 

molecules 

uv, visible 

2-20 










Eaman scattering cross sections. This physical fact leads to the immediate 
conclusion that the most sensitive remote detection methods are hy absorption. 
Resonance fluorescence and fluorescence also have cross sections that are 
nearly that of absorption. However, for monitoring in the troposphere 
quenching reduces the fluorescence cross section. Quenching is not a problem 
at high altitudes and remote monitoring by resonance fluorescence has a clear 
advantage in that case. 

We can now return to Eq, (8) and consider the effective reflectivity 
of the "target” for remote sensing. Table III lists the scattering process 
the corresponding effective reflectivity. 

For absorption measurements with a remote sensor or retroreflector and 
a collimated transmitted beam, the range squared dependence is effectively 
cancelled. Thus all of the transmitted power is collected except for a 
collection efficiency factor g , which is near unity, and the effective 
target reflectivity is orders of magnitude better than for a Lambertian 
scatterer. Letting p/jt = R |/A in Eq. (8) reduces it to the form commonly 
used to describe long path absorption. , 

For long path absorption by scattering from topographic targets, p 1 
for aLl infrared wavelengths and p .1 for ultraviolet and visible wave- 
lengths. Thus the effective target reflectivity is approximately l/jt . For 
single ended long path absorption measurements from topographic targets over 

4 6 

kilometer ranges, 10 - 10 more transmitted power is required than for doubled 

ended measurements using a remote detector or retroreflector. Since the required 
power for double ended long path absorption is less than one microwatt,^ single 
ended measurements can be made wit*h milliwatt to watt laser sources. The 
HeNe laser geodolite ranging systems are an example of single ended range 
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TABLE III 


EFFECTIVE EEFLECTIVITY 


Method 

(pA) 

Magnitude 

Long- Path Absorption 
with Retroreflector 


10*^ - 10^ 


A 


Long Path Absorption 

p/tc 

lA 

from Topographic Target 



Differential Absorption 

/n\ PMie 

10 ^ 


\2 / 


Resonance Fluorescence 
and Fluorescence 

^ 

10 ^ (l ppm) 


\2/ 

(uv and visible) 

Raman 

/cr\ Ncr 

/ j RAM 

10 (l ppm) 


\2/ kn 

’ 10‘^ (1 atm) 
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measuresients by scattering from topographic targets. 

Remote detection by backscattering involves the range resolution ~ 
as well as the backscattering cross section • For Kie backscattering the 
returned power is independent of the pollutant concentration and depends only 
on the Mie backscatter coefficient 3,,. . For resonance fluorescence and 

Raman backscattering, the returned signal depends on the scattering cross 
section and on the pollutant density. The dependence on pollutant density 
reduces the returned signal for weak pollutant concentrations and thus also 
reduces the measurement S/n ratio which further effects the detection 
sensitivity. For depth resolved pollutant measurements differential absorption 
has the advantages of a relatively large effective reflectivity which is 
independent of pollutant concentration and the additional advantage of the 
large absorption cross section which is not affected by quenching. 

To complete the general formulation of remote pollutant detection we 
need to introduce detector parameters and evaluate the detector S/N ratio for 
a given returned power. This procedure has been discussed previously, 
but a summary is useful here. 

The voltage signal to noise ratio, for a detector is 

where is the received signal power and if is the amplifier bandwidth. 

The first term in the denominator is due to the shot noise generated by the 
signal itself. Here Iw is the photon energy and is the detector quantum 
efficiency. The second term in the denominator represents the noise of the 
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detector in the absence of the signal (dark current and background noise and 

preamplifier noise characterized by the noise equivalent power NEP). A factor 

of 2, added under the square root, takes account of the noise resulting from 

two measurements required to determine the pollutant concentration. This noise 
2 2 2 

is s (p.,) + (P„) where on and off refer to making the measurement on 

N N on N off . . ® 

resonance and off resonance of the characteristic line of the pollutant. Usually 
either the first or the second term in the demoninator need to be considered. 

For shot noise limited detection the dark current term in the denominator 
can be neglected. For practical purposes this limit can be attained with photo- 
multipliers if all background radiation is adequately filtered. At a given 
signal to noise ratio, the minimum detectable power from Eq. 10 is given by 

I^^“(Ehot noise) ^ l^(^) r (|)^ if (U) 

where the factor F has been added (between 2 and 5) to account for the partition 
noise effect of the multiplier dynodes. Equation (11) gives the minimum detectable 
power for current measurements. 

Riotomultipliers can also be used in the photon counting mode. For count 

3 “1 

rates greater than approximately 10"^ sec the photon counting method and current 
detection method give equal sensitivities. However, for low count rates photon 
counting offers an improvement in signal to noise ratio by the factor F . Thus 
for weak signals photon counting is preferred while for more intense signals and 
vide dynamic range the boxcar integration technique is more useful.^ 

An important extension of the shot noise detection regime is heterodyne 
detection. The nonlinear response of the photo- detector allows the generation 
of a mixing frequency at the intermodulation frequency of the two applied light 
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fields. The heterodyne signal to noise ratio for detection at the IF frequency 


is 


62 


(heterodyne) > 


2 hv(|) Af 


(photoemitter reversed (l2a ) 
biased diode) 



Af 


Tl 


( phot oc onduct or 
photovoltaic 
diode) 


(12b) 


where the extra factor of 2 for the photoconductor detector is due to 
generation-recombination noise. 

The required local oscillator power is such that the shot noise term in 
the denominator of Eq. (lO) equals the dark current or amplifier KEP term. 

This condition can be written in a number of ways depending on the expressions 
for -the HEP. However, if amplifier noise dominates and is expressed as an 
equivalent noise temperature T. then 

A 



(k T^) hw 

2~2 
T| e G 


( 13 ) 


where e is the electron charge, G the detector gain and R,. the load 

L 

'■ 

resistance. For example, for an amplifier with noise temperature of 2k0 K 
(noise figure 2.h dB) and assuming = 0.5, G = 0.12, = 50 n and \ = 10 uirf 

we obtain = 12 mW . This is a rather high local oscillator power since it 

may result in detection heating. However, wide band photodiode detectors with 
\mlt quantum efficiency should require much less (~100 m W) local oscillator 
power to reach the quantum limit. 
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Heterodyne detection offers extreme spatial filtering due to its antenna 

like properties and the short optical wavelength involved. Siegman has 

. ^ ^ 2 

shown that the effective aperture of a heterodyne detector is ^ 

where Aj^ is the effective receiver aperture and 0^ the receiver solid angle. 
Thus in addition to orders of magnitude more sensitivity than darh current 
limited detection^ the heterodyne detector offers diffraction limited spatial 
filtering against the incident background radiation. 

A potentially useful form of heterodyne detection for remote monitoring 
is infrared up-conversion in a nonlinear crystal. The noise power for an up- 
efficiency F onto a shot noise limited photomultiplier is 

- 2 

Af 

‘ (14) 


4 hv(:g) Af 
(up-conversion) = 




which is equivalent to the previous expression for heterodyne detection by a 

photoconductor. Using a cw argon ion laser at IW power to pump a LillbO^ 

68 

crystal Smith and Mehr were able to up-convert and detect 3 ^im radiation at 
a HEP of 10 W . Up-conversion offers wider acceptance bandwidth (*^100 GHz) 
comj^red to heterodyne mixing detection (-^IGHz) and yet still maintains spatial 
and spectral filtering. 

The useful properties of heterodyne detection has led to the suggestion of _ 
remote detection of pollutants by monitoring thermal induced fluorescence from 
hot effluent gases 


Menzies^^^"^^ has evaluated the use of heterodyne detectic 


•51 

for remote monitoring and has recently-^ demonstrated room temperature detection 
of ozone, sulfer dioxide, ammonia and ethylene using a CO^ laser as a local 

oscillator. 
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Heterodyne techniques can also be applied to improve the detection 
capabilities of any remote monitoring method. However, the cost and complications 
of present tunable coherent sources has prevented the wide application of hetero- 
dyne detection. 

For dark current or background limited detection the minimum detectaole 
power evaluated from Eq. (lO) is 

I^°(dark current) = HEP (^) \/2if , (15) 

This type of detector noise limitation applies, in particular, to photo- 
conductors and other detectors used in the Infrared. It is characteristic 
for dark ciirrent limited detectors that they achieve, for constant average 

pover, a higher signal-to-nolse ratio with increasing peak power in pulsed 

71 

operation. 

* 

Detector performance is also expressed in terms of the detectivity D 
which is related to the HEP by D* = ■sA/HEP where A is the detector area. 

* 10 "i- 

Typical values of D are 6 X 10 cm Hz-^/W for room temperature silicon 
diodes, 3 X 10^ cm Hz^/W for 77°K operation of InAs (1 pm - 3.3 pm) photo- 
voltaic detectors, 1.1 X 10^ cm mP /M for 77°K InSb (I .5 - 5-5 m) 

photovoltaic detectors, to 2 X 10^^ cm for 77 K HgCdTe(2 pm to 25 pm). ~ 

detectors. The rise time of these detectors is less .than 100 nsec. Finally, 
detector costs vary from approximately $100 for silicon diodes to $1,000 for 
recently developed HgCdTe detectors. A complete manufacturers reference to 

72 

detectors is given in a recent trade Journal issue. 

In the above expressions fof detector minimum detectable power, we have 
assumed Af is the electronic bandwidth required to resolve a pulse. More 
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specifically we take tfr =2 as the condition necessary to resolve a pulse 
of width T • For pulsed generation the effective detector bandwidth decreases 
with the number of averaged pulses n as Af/n . 

We are now in a position to combine E-q,- ( 8 ) for the received power with 
the expressions for minimum detectable power and solve for the detection para- 
meter of interest at specified values of the other variables. The parameters 
of interest are required transmitted energy vs range, s/u ratio vs range 
at a fixed energy, and minimum detectable pollutant concentration relative to 
one atmosphere , t] ^ at a given transmitted energy and range. The above procedure 
is carried out and the results are discussed for Raman scattering, fluorescence 
backscattering, the integrated path absorption methods, and for the differential 
absorption method in the following sections. 

\ 

' B. Raman Method 

Remote atmospheric monitoring by Raman backscattering was first demonstrated 
20 

in 196’J by Leonard. The early work in Raman scattering including the detection 
21 22 2*5 2i^ 

of water vapor * ^ and SOg and COg was discussed in the introduction. 

Early remote pollutant detection work was limited to Raman scattering 
because of the lack of high peak power tunable laser sources. The principle 
advantage of the Raman method is the use of a fixed frequency laser transmitter. 

The other apparent advantages include operation in th^ visible or ultraviolet 
and the use of photomultiplier detectors and the detection of homo-nuclear 
molecules addition to other molecules of interest. On the other 

hand, the disadvantages of remote detection by Raman scattering include very 
small Raman scattering cross sections and therefore lack of sensitivity, operation 
in the visible with high energy laser sources and the resultant eye safety problem. 



and simultaneous excitation of the Raman scattering of all molecules present 
thus requiring a spectrometer filter to adequately isolate the Raman signal 
of interest. Of these disadvantages, it is the lack of sensitivity that has 
forced thinking about improvements of the Raman method and the remote detection 
of pollutants by other methods. 

An interesting benefit of the remote Raman detection work has been the 

4 

measurement of absolute R^man cross sections of gases. Inaba and Kobaysi 
provide a table of measured Raman cross sections in their review paper. 

7*3 

However, additional measurements are still being made as evidenced by a very 

74 

recent absolute Raman cross section measurement of 

To evaluate the sensitivity and required transmitted energy for the Raman 

method we substitute the expression for the minimum detectable power for a shot 

noise limited photomultiplier given by Eq. (11) into Eq. (8) and replace (o/at) 

Kij / \ ^ ^ 

in Eq. (8) by ■ ■ ( '!”) given in Table III. Solving for the transmitted 


powear we find 


hit 


P - 

0 K 


4(^) F (|) Af 
T} N 


4it 


lij 


RAM 


2 R 

— — exp 
TC A 


R 


/ 


dr 


(16) 


where the terms in the brackets give the detector minimum detectable power 
for a shot noise limited detector. For shot, noise limited detection the highest 
average power is of interest. However, for depth resolution and to remain in 
the shot noise limited detection range, pulsed laser sources are usually used. 
Multiplying Eq. (l6) by t , the required transmitted energy of pulse width t 
is given by 
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As an example, ve consider the Raman detection of CO at 0.3^7 jjm, the 

doubled Ruby laser output wavelength. At this wavelength (dc/dO)^Q = 

““^0 2 ^ 

3-6 X 10 cm . We let the atmospheric volume extinction coefficient be a 

constant 1 lun ^ corresponding to a visibility of approximately 5 leu, and assume 

a photomultiplier detector with F = 4 and a 20 MHz bandwidth required to 

resolve a 100 nsec pulse. For a U , 100 nsec transmitted pulse and a 
2 

1000 cm area receiving telescope at K = .1 we find the minimum detectable 
pollutant concentration vs range shown in Fig, 2- Also shown in Fig. 2 are the 
peak and average CO levels for urban atmospheres. 


It is readily apparent that the Raman method does not have the sensitivity 
required for ambient pollutant concentration measurements. However, the Raman 
method does have adequate sensitivity to monitor , Og and 1^0. to a few 
kilometer range. 

Using a nitrogen laser operating at 3371 ? at 20 kW peak power 10 ns half- 

7'5 

Width at 50 pps, Kobayasi and Inaba'"^ have detected Raman signals from clear 
air and from an oil smoke plume. The receiving system consisted of a 30 cm 
diameter telescope followed by an f 8.5 half meter grating monochromator and spectre 
interference filter. The signal was detected by a photomultiplier and boxcar 
integrator sampling circuitry and displayed on an X-Y recorder. Figure 3a 
shows the return from ordinary air at 30 m range. The Ng and Og, 0 and S ~ 

as well as Q branches are evident. Also evident is the strong Rayleigh and 
Mie backscattered signal which must be suppressed for Raman measurements but 


which is very useful in the differential absorption method considered later. 

Figure 3b illustrates the Raman return from oil smoke. Here the Ng and 
Og peaks provide a calibration for the other molecular return signals of interest. 



However, the spectral width and intensity of the 0 and S branches tend to 
interfer with the Raman signal of the more wealcly concentrated smoke plume 
constituents. Also indicated is the fluorescence return which may be a useful 
signal for some remote monitoring applications discussed in the following 
section. 

The use of a photomultiplier detector, and the intense Rayleigh and Mie 
backscattered signal compared to the Raman signal presents significant filtering 
problems against background signal. The requirements for background filtering 
were discussed by Kildal and Byer who showed that the background signal must be 
less than the Raman signal if sensitivity is to be maintained. This requirement 
forces the consideration of grating monochromators and interference spike filters 
in the receiving system thus adding complications and lowering its optical 
efficiency. Recently, Cooney proposed a differential Raman measurement using 

Ng and Og as a reference to eliminate the common background and allow. daytime 
operation. 

progress in laser sources has led to improvement in Raman detection systems. 

77 

For example Hirschfeld et al. , ' describe a frequency doubled Ruby laser trans- 
mitter with a 3 m receiving aperture telescope system which has been used in 


the field to detect COg , 1^0 and SOg over a 200 m range. The 2J, 2pps Ruby 
laser source was doubled in KDP with 8^ efficiency and the returned signal was 
processed by photon counting techniques referenced against the nitrogen signal 


in a differential manner. Range gating provided a 10 m depth resolution. This 
Raman detection system probably represents the state of the art for sensitivity 

78 

end range. 


Recently Raman scattering has been used to vertically probe the atmosphere 
to a height of up to 6 km. The Raman-laser system is described by Kent et al. 
and the measurements and probably extension to hO km height is discussed by 
Garvey end Kent.®*^ Xhe application of Raman scattering return from for 
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temperature measurements has been considered by Strauch et al. Experimental 

results over a 30,5 m range at a 5 m depth resolution confirmed the R-=man ^ck- 

scatter vs temperature correlation and indicated a differential temperature 

o 

sensitivity near 3 . 

Despite the above success of the Raman method^ its lack of sensitivity 

has severely limited its application to ambient pollutant monitoring. A 

number of suggestions have been made to improve the Raman sensitivity. They 

include using resonance to enhance the Raman cross section, using interferometric 

methods to improve receiver entendue and applying four wave Raman mixing to 

generate a coherent anti-Stokes signal vd.th high efficiency. 

Resonance Raman scattering has achieved considerable recent attention 

due to the potential of improving the Raman cross section by tuning the probing 

laser radiation near an electronic resonance. The general form of the third 

X gT 

order susceptibility in Eq. ( 5 ) is° 

^( 3 ) « ^ 

K - + i - CO;,) 

where cd^ and o:^ are electronic frequencies and is the Raman frequency. 

As the exciting frequency ^approaches an electronic resonance 

is resonantly enhanced. The enhanced Raman cross section still falls well shorty 

1 82 

of the on resonance fluorescence cross section. ^ Experimental work has 

83 8^-8 t 

verified the theory for both atoms ^ and molecules. ° For remote detection 

the tunable laser source required to achieve Raman resonance enhancement is 
probably better utilized by tuning exactly on resonance and exciting resonance 
fluorescence. In that case the fluorescence return, even if quenched, is still 
about two orders of magnitude greater than the resonance Raman return.^ 
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The . rotational Raman scattering cross sections are larger than the 

■vibrational-rotational cross sections and interferometric filtering has been 

88 

proposed to take advantage of this to improve Raman sensiti-vities. However, 

the rotational lines of most molec\iles lie within 100 cm ^ of the exciting line 

and overlap interference is a severe problem. These reasons have prevented 

8q 

realization of the expected improvement in sensitivity, ^ 

A novel application of tunable laser sources to Raman detection is 

At 

coherent anti- Stokes Raman spectroscopy by four wave mixing. In this 
coherent Raman scattering process, a tunable laser mixes with a pump laser 
to generate a coherent laser-like anti- Stokes output beam. The Raman suscept- 
ibility shows resonance when - o:^ - = 0 according to Eq. (l8) where 

is the pump frequency, the tunable frequency and is the Raman 

resonant frequency of the medium. 

The coherent antl-Stokes beam is generated with a conversion efficiency 
given by-^-'^ 

P(a) ) 2.77 X 10“^ 

^ 2 1^ ^rAM 

irtiere N is the density, is given by Eq. (5) in terms of 

^coh the coherence length in cm which in liquids is approximately 
0.1 cm but in gases is near 100 cm, and is the pump intensity in W/cm^. 

Here N is the bulk Raman susceptibility in esu units (cm^/erg) and 

is the anti- Stokes wavelength in cm. 

The advantages of four wave Raman mixing are the very high conversion 
efficiencies to the anti-Stokes frequency, the spatially coherent output beam, 
and the discrimination against background fluorescence. These advantages have 



■^coh ^1 


(19) 



been experimentally verified in liquids where up to 1^ conversion efficiencies ‘ 
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have been observed in benzene. 

Four t/ave Raman mixing has also been used to measure small gas concentrations 
91 

by Regnier and Taran. They used a Ruby laser pumped stimulated molecular 
hydrogen Raman source mixing with the remaining Ruby light to measure 
concentrations down to 100 ppm in nitrogen. Below 100 ppm the background 
third order susceptibility of gave rise to a signal. which interfered with 

the weak anti-Stokes signal. With presently available high peak power 

tunable laser sources four wave Raman mixing should provide a valuable 
technique for the investigation of gases down to the ppm level. Although 
the input beams must be focused in the medium, the spatially coherent output 
beam can be remotely monitored thus providing a local pollutant monitoring 
system with a probe that is completely free from the volume of interest. Recent 
work has verified the important parameters of four wave Raman mixing in liquids 
and gases, but more research is necessary to evaluate the remote detection 
possibilities. 

C- Resonance and Fluorescence Scattering 

92 

In 1969 Bowman et al.,^ used a tuned dye laser to detect sodium in the 

upper atmosphere by resonance backscattering. By 1970 regular measurements^^ 

were being made of the sodium double layer at 90 altitude. Since that time 

95-98 

observations of the sodium layer have been numerous and recently they have 

99 • loo 

been extended to daytime^ and to potassium. 

The use of resonance and fluorescence backscattering for probing the 

troposphere has been limited. However, the method has been extensively analyzed 

by Kildal and Byer^ and recently by C.M. Penney et al.,^^^ Resonance backscattering 

75 5 

has also been discussed by Kobaysi and Inaba and by Measures and Pilon."^ 



To evaluate the returned signal power due to resonance and fluorescence 
hackscattering we proceed as before for the Raman method. However^ we must 
now consider the additional complication of a finite fluorescence lifetime 
which both reduces the returned power and smears out the depth resolution. 

An exact analysis of the returned power for a finite fluorescence lifetime 
was given by Kildal and Byer^ and is summarized here. 

In general, the returned signal power is given by 


P^(E) = KE^ — 

sp 


,«r 


s(t) 


^min 


( 20 ) 


where oOp and are the fluorescence and transmitted frequencies, is 
the fluorescence cross section without quenching , is the range to the 

beginning of the pollutant cloud. S(t) is a dimensionless integral with a 
peak value of unity given by 


T 

P 


?"f/ 2 dE / 

1 SR\ 

/ - 2(1 -exp 

- - ft - -1 

•'dt-T )/2 ^ \ 

L ^ ^ '1 


• exp[-aj,N(R - R^.^)] + [l - expC-T^T^)] J 


=i*-"pV2 dE 
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where is the radiation decay time which is less than t , the 

F sp 

spontaneous decay time^ due to quenching and is the laser pulse width. 

The integral S(t) has been studied in detail and it is plotted vs time 

for a number of cases of interest,^ In general it reflects the increased 

time over which fluorescence returns to the detector due to the finite 

radiation decay time of the fluorescence. The fluorescence decay time leads 

to a depth resolution given by Eq. (9) in which t is the dominant term. 

r 

For infrared transitions the depth resolution varies between .15 km and 15 km 
thus preventing depth resolved measurements. For visible and ultraviolet 
transitions^ quenching reduces the radiative lifetime to the order of nano- 
seconds so that the depth resolution equals that for Raman scattering. 

For visible and ultraviolet transitions, the integral S(t) simplifies 


to*. 


= f ) - 


CT„ exp[-a^N(R - 


2 K(R - c T^ 2 ) 


If we let R » c T ^2 which is the usual case of interest, then S(t) 
reduces further to 

s(StravLlet) “ \in ^ ^ 

where we have incorporated the pump depletion due to absorption by the fluorescing 
species exp -o.^(r - R ) into the total atmospheric extinction term 


exp 




H 

f a (r)dr 
0 ^ 


] 
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Substituting Eq. (22) into Eq. (20) gives for the power received at the 
detector due to backscattered, fluorescence 


P^(R) = KE^ 




c ^ 
2 C, 




f q: (r)dr 
0 ^ 


. (23) 


If we note that a 


ABS 


= 


where Q is the Stern- Volmer quenching factor 
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with l/Z the probability of quenching per collision and t , the collision 

col 

time which is approximately 10 ^ sec at atmospheric pressure^ then Eq. (23) 

canvbe written in a form identical to that for Raman scattering except that 

Q . Thus . ' ' 


^RAM replaced by a 




ABS 


KP 


ABS 


o 


4jt 


cT cv A r ^ 1 

Qf— — ^ exp -2 /a (r)dr 

2 0^ L 0 ^ J 


(24) 


where the effective reflectivity 



equals 


scattering in agreement with Table III, 



for quenched back- 


The two cases of interest in evaluating fluorescence backscatter are 
infrared backscatter from molecular vibrational- rotational transitions and 
visible or ultraviolet resonance and fluorescence backscatter from atoms and 
electronic transitions of molecules. Fluorescence cross sections and quenching 
factors have been listed and are not repeated here.^^"^^ 
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For. the infrared case the power received at the detector given by Eq.; (20)_ 
is equated to the minimuni detectable power for a dark current limited detector 
given by Eq. (15) and the minimum pollution concentration is found as a function 
of transmitted energy and range at a given set of parameters. In this case, 
only a fraction, f , of the total fluorescence band is monitored, for example, 
a single vibrational- rotational line. The integral S(t) cannot be simplified 
and must be evaluated. Figure 4 shows the result of an evaluation for CO at 
4.6 |im assuming an InSb detector, a 100 mj transmitted pulse and a S/N = 1 . 

The other parameters of interest are given in the figure. The integral S(t) 
evaluates differently for thin and thick pollutant layers. For a thin layer 
the pump depletion due to absorption by the pollutant is negligible while for 
a thick layer the transmitted pump beam is totally depleted within, the layer. 
Optically thick pollutant layers prevent depth resolved measurements due to 
transmitted pump depletion over a length large 

atomic and molecular electronic cross sections make optically thick pollutant 
clouds a distinct possibility in the visible and ultraviolet. Measures, and 
Pilon^ present calculations showing the resonance detection of optically thick 
layers and the resultant pump depletion and loss of depth resolution. 

For the visible and ultraviolet case the power backscattered to the 
detector is given by Eq. (24). Equating the received power to the shot noise 
limited minimum detectable power given by Eq, (ll) gives for the required 
transmitted energy 

where f is the fraction of total fluorescence monitored, = ct /2 is the 
range resolution and the other quantities have been previously defined. 
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Equation (25) is similar to the previously derived result for Raman scattering 

given by Eq. (1?), except that the required transmitted energy is reduced . 

■ -2 

from the Raman case by the ratio which is approximately 10 

“8 —10 

for molecular pollutants and 10 - 10 . for atomic pollutants such as 

sodium and mercury. Examples of detection sensitivities for electronic 

molecular transitions and atomic transitions are given hy Kildal and Byer, 

7®5 . 5 

Kobaysi and Inaba and by Measures and Pilon. 

Recently Gelbachs et al. utilized resonance fluorescence of NOg 

excited by an argon ion laser to detect NOg , locally with a sensitivity of 

one part per billion. The laser excitation was at .488 pm and the fluorescence 

was monitored at 0.7 to 0.8 pm . Local Los Angeles air, drawn through filters, 

was monitored and HOg variations measured over a few hour period showed peak 

concentrations near 0.1 ppm on smoggy days and O .03 ppoi on smog free days. 

Filtering the air was necessary to eliminate interfering fluorescence from 

particulates in the air. 

The Increased sensitivity of resonance and fluorescence backscattering 
makes it appear useful as a remote monitoring method compared to the Raman 
method. However, the requirement for a tunable laser, the lack of depth 
resolution in the infrared, and the not precisely known fluorescence cross 
section due to the unknown local quenching factor, severely limit the usefullness 
of resonant backscatter or fluorescent backscatter as a remote pollutant measure- 
neht method in the troposphere. Fina31y, as will be shown later, the sensitivity 
of the fluorescent backscatter method is surpassed by that of the differential 
absorption method which does not suffer from quenching or return signal strength 
dependence on pollutant density or cross section. 
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An alternate approach to monitoring radiated fluorescence is to apply 

the tunable coherent radiation to the detector to increase its sensitivity by 

heterodyne detection. Heterodyne detection was discussed in Section IIA so 

here only an estimate of detection sensitivity is considered. 

‘ -2 -1 -2 

The sky brightness in watts m Hz rad due to thermal radiation with 
emissxvity for a single polarization is ■ 

B=(hv3/o^)( ) , ( 26 ) 

\ exp tw - 1 / 


The power received at the heterodyne detector is given by Eq. (?) with 
AX replaced by the I.F. bandwidth^ B , and ^ heterodyne 

detection^^ so that 


% 



exp(b\)/klg) - 1 


( 27 ) 


vhere is the background atmospheric emissivity. 

B 

The signal received from a pollutant cloud filling the field of view 
at temperature is 


[1 - exp(-Kj,^„ri)] 

^ IwB 

^ exp (i^j/kT^) - 1 


( 28 ) 


where the emissivity of -the pollutant gas is related to the absorptivity by 
Kirchoff's law. The signal to noise for heterodyne detection of an emitting 
pollutant cloud is found by setting (P°^ - P°^^) equal to the heterodyne 
minimum detectable power given by Eq. (I2b). The resulting voltage signal to 
noise ratio squared is 


2 

\ H 

1 - exp(-Na^ggL) 


I " k 

exp (hv/kTp) - 1 

exp (hv/kTg) - 1 


Bt 


( 29 ) 
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where r = 1/Af is the integration time. Since the emissivity of the 

background is usually less than 1, the temperature of the pollutant gas 

can be less than the surrounding air and still be detected. 

. Figure 5a aiid 5b show the sensitivity limits of heterodyne detection 

of SOg and COg vs the IF bandwidth- integration time product. The 

detectable concentration is given as an equivalent absorption coefficient 

where 10 atm cm is equivalent to 1 ppm km. Figure 5a shows the sensitivity 

70 

limits calculated by Menzies for SOg with a measured absorption coefficient 

“1 “1 T p jQ 

of 0.58 atm cm at room temperature using a C Og laser as the local 
oscillator. Figure 5b. shows the calculated COg detection sensitivity at 
T = 450°K for all the mixing signal falling on the detector (solid line) 
end for 2C^ falling on the detector which corresponded to the experimental 
condition. The crosses indicate actual measured sensitivities. 

These results show that heterodyne detection is useful for remote 
monitoring of hot effluent gases from concentrated sources. However, the 
sensitivity of the heterodyne method is still not adequate for dispersed 
pollutant measurement and is not nearly as sensitive as absorption methods 
to be discussed in the following section. 

D. Long Path Absorption Methods 

Long path absorption for the detection of atmospheric pollutants using 
laser sources is an extension of well known methods. However, the 

laser source with its spatial coUimation, very high spectral resolution and 
high power offer significant advantages over previously used incoherent sources. 
Closed path double ended long path absorption and open path single ended long 
path absorption using topographic&l targets are discussed in this section. 

These measurement methods provide the most sensitivity at the least required 
transmitted power of any of the remote detection methods. 
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Bemote air pollutant detection using laser sources was discussed as 

106 28 107 

early as I967 by Jacobs and Snowman. Later Hanst ' described pollutant 

detection using fixed frequency laser coincidences with absorbing lines. More 

recently, Zaromb,^^ and Japanese worker described remote sensing by 

laser absorption and Hodgeson et al. considered advanced monitoring 

techniques. Kildal and Byer^ have given a detailed analysis of the doubled 

ended long path absorption method. The analysis was later extended to the 

single ended absorption method using topographical targets by Byer and 
33 

Garbuny. 

Double ended long path absorption experiments have been carried out 

using low power diode lasers by Hinkley^^ ^ end using COg lasers by Snowman. ^ 

Ilk 

Recently an argon ion laser has been used to detect NOg. The topographical 

11 ^ 

single ended absorption method was first demonstrated by Henningsen et al. , 
in the. infrared using a parametric tunable laser source. The parameters of 
Importance for long path absorption measurements and the detection sensitivity 
are discussed in this section along with recent experimental results. 

The long path resonance absorption method measures the total integrated __ 
pollutant concentration over the path. It has the advantages of increasing 
sensitivity with range, the largest interaction cross section which is not 
modified by quenching, and the least required transmitted laser power. It 
has the disadvantages of being double ended and of lacking depth resolution. 


/ 
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The transmitted intensity at the detector follows from Eq. (8) if- we 
let (p/jt) = R^/A , assume a round trip, and separate the beam attenuation due 
to pollutant absorption from the atmospheric attenuation. We then have 

P°"(R) = KP^exp^-2/ a^sH(E') dH'j 

. exp |^- 2 ra^(R') dE'J (30) 

vhere N(R' ) is the pollutant density over the range R and we assume that 
the transmitted intensity is less than the absorption saturation intensity. 

The returned power for the pump frequency tuned off the absorbing 
transition is 


V 





yor-CR ) dR' 
0 ^ 


] 


(31) 


so that the integrated pollutant concentration is determined by the log 
ratio of the on and off returned power as 


R 

2 / N(R')<iR' 
0 



(32) 


Assuming that it is possible to measure a 1^ change in transmitted power 
when tuned on and off the absorbing line, the minimum integrated density is 
given by 



/N(R') dr 

0 J min 


0.01 


^ABS 


( 33 ) 
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Table IV lists representative absorption cross sections for molecules 
and atoms and gives the minimum detectable pollutant concentration for an 
absorption length R of 100 meters. Since the absorption cross sections 
are the largest of the Raman, resonance or fluorescence cross sections 
these sensitivity limits are the best obtainable by optical means. Thus 
there is a physical limitation for pollutant detection that leads to a depth 
resolution- sensitivity trade off. Fortunately, ^ABS lai'ge enough that an 
absorption distance of 100 m allows ambient pollutant levels to be measured 
to less than 20 ppb. - 

The absorption sensitivity can be improved by measuring less than a 1^ 

signal power change. The allowed measurement accuracy is discussed by Kildal 
1 

and Byer and by Byer and Garbuny.*^"^ These authors show that to detect a 

-7 

0.1^ change in intensity with an InSb detector requires only 10 watts of 
incident power at the detector. There is an optimum value of absorption cross 
section given by N = 2.22 and IT = 1.3.1 for shot noise and dark 

current limited detection. In addition, the signal to noise ratio required 
to achieve a measurement accuracy Ax/x is 

S/K = Ax/x (3k) 

where x = 2 ■ 

Other parameters which affect the measurement accuracy and limit the 
detection sensitivity are atmospheric turbulence and laser power fluctuations. 

For Ctf measurements the atmospheric turbulence is a significant limitation 
for chopping rates on an off resonance of less than 200 cps. Since the signal to 
noise ratio improves with higher peak power for dark current limited detectors, the 
use of Q-switched or pulsed laser sources is advantageous. The laser power 
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TABLE; IV 

Minimum Measurable Concentration for an Absorption Length of 100 Meters 



CO 

NO2 

S O2 

C6«6 

Na 

Hg 


4.7 

.4 

.29 

. 29 

.5896 

. 2537 


l.8x 10''^ 

0^ 

1 

0 

X 

00 

CM 

3. 4x10" 

1.3x10"'^ 

4.8x10"'^ 

5.6x 10"'^ 

77 ( 1 %) 

2. 1 X 10"^ 

1.3x10"^ 

1. 1 X 10"^ 

1. 5 X 10"^ 

7.9x10"'"^ 

6.7x10"'^ 



^actuations can be reduced by careful engineering or by ratio normalization 

of each laser pulse. The ultimate accuracy achieved by these steps is 

determined by the dynamic range and linearity of the detector. However 

0.1^ ratio accuracy is routinely achieved. 

The required transmitted power at a given signal to noise rat-io for 

33 

measurement of a pollutant concentration has been shown to be 





(l)= 

>(rf7p,)' 

(shot noise) 

(l)= 


(dark current) 

s 

N “ 


(shot noise) 

S 

N 

•35te"/P,) 

(dark current) 


(35) 


where is given by Eq. (11) for the shot noise limited detection and by 

Eq. ( 15 ) for the dark current limited detector with S/N * 1. 

For molecules x varies between .1 and 1 for R = 100 meters at a 
concentration of 1 ppm. For atoms x varies between 10 and 3 OO 
R ss 100 meters at a concentration of 1 ppm. Assuming x > 1^ for the detection 
of CO at i}-.7 }im with an InSb detector, the required transmitted power from 
Eqs. (31) and (35) is 




(36) 
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For an MEP = ,3 x 10 ^ , K = .1 a bandwidth of 1 MHz , ^ = 100 

which implies a detection sensitivity, and assuming ~ 1 km ^ ^ the 

required peak transmitted power for a 1 km range is = 1 mW in a 

2 psec'pulse. This is well within the range of cw and pulsed infrared lasers 
112 

and diode lasers. 

We can invert the above procedure and ask what is the maximum range for 
a 1 kW peak power laser source. For = 1 km ^ we find R = I 6 km at 

g 

a (— ) = 100 , Beyond that range the signal to noise ratio decreases reaching 
unity at 20 km. 

These examples show the advantages of long path absorption Measurements 
which include low required transmitted power, long ranges and of course good 
sensitivity. However, in some situations the requirement for a double ended 
system is not possible to meet. In those cases the single ended topographical 
absorption technique is an alternative measurement method. 

However^ if the .transmitter is tunable, then it is possible to sample the optical 

path for absorption from specific pollutant molecules. The measured integrated 

density is again given by Eq. (33) ^ measured intensity change. ' Table 

III again gives the detection sensitivities for a 100 m path. 

The advantage of absorption measurements using topographical targets is 

the single ended measurement with the sensitivity of the closed path absorption 

method. The disadvantage is the requirement for increased transmitted power. 

To find the required transmitted energy for absorption using topographical 

targets we must consider the measurement procedure carefully since pump depletion 

33 

due to the pollutant itself may be important. 
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The power received off resonance is given by Eq.. (0) with an effective 
reflectivity (p/s) , thus 


■ [(;) ” 4 ] 


exp - (2 R) 


(37) 


iriiere cc » or -fc a is the atmospheric extinction coefficient given by 
Eq. (2) without the absorption term. If we assume that the atmosphere may 
itself contain the pollutant gas in a normal concentration which is less 

than the pollutant concentration , then the power received on resonance is 


given by 


[“so ^ "'a °ABS 


exp 1-2 a 


ABS 


{ v» 


') dH'j 


(38) 


where L is the extent of the pollutant cloud. The difference in the power 
received off resonance and on resonance found by combining Eq. (37) and Eq. ( 38 ) 

is 


r r r 


1 - exp [ -2 + IT L ] o 


ABS 


(39) 


This difference must be measured relative to the power received at the detector 
but attenuated by the normally present concentration of gas in the atmosphere 
given by 


P - P 
r r 


normal 
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1 - exp |-2 N^R 
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For a measureaient of the pollutant concentration the difference between 
(39) Eq. (4o) must equal the noise power given by Eq. (I5) at a 
given signal to noise ratio, or 

Pp exp[-2 N^R | ^ " axp[-2 N^L | = ^ 



( 41 ) 


Substituting for 
power gives 



from Eq. (37) and solving for the required transmitted 



exp -f g^g] R 

1 - exp (-2 NpL 


1 

Finally multiplying the above result by t and setting Afr » 2 yields for 
the required transmitted pulse energy 



I H (1)'^]? 


f ^[”se ^ 

A 1 - exp(-2 K^L a^g) 


(te) 


•aq 

a result first derived by Byer and Garbuny. The required transmitted energy 
varies as as expected for a dark current limit detector. It also increases 

with range due to the Lambertian scattering properties of - the topographical 
reflector with reflectivity p into 2-r stearadians. 

Equation (42) gives the required transmitted energy assuming that the 
pollutant of interest is also normally present in the atmosphere. The equation 
simplifies considerably for most pollutants which are normally not present in 
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the atmosphere. For this case, the term exp(2 R) approaches unity. 

Finally, if the pollutant cloud itself is optically thin, then the denominator 
6 f Eq. (42) reduces to 2N^L . In the limit of a very tenuous pollutant 

Cloud the reduced value of 2N L leads to a higher required transmitted 

P iU30 

energy. Thus we expect that there must be an optimum absorption cross section 
for a given pollutant density and cloud thickness, 

Byer and Garbuny have investigated the optimization of c for the minimum 
required transmitted energy. Setting dE^/da =s 0 In Eq. (42) leads to 

'’min = ^ V 

Where q = N L/lf,R . If N L < < N.R , a situation that may hold for CO 
p ' A P 

detection, then o , = 1/2 N.R . In other words at the optimum value of a 

' min ‘ A 

2R equals the e"^ depletion length. For example, CO normally present 
in the atmosphere at 0.2 ppm concentration, requires at exact resonance ■ 
a s o^Av ^ transmitted laser energy 10 ^ larger than for a = 
measurement over a 5 range. 

As an example of long path absorption measurement we consider the detection 
of CO at 4-. 7 pm using both long path absorption methods. Consider first the 
required transmitted energy for topographical target scattering given by Eq. (42). 
We assume an InSb detection with an HEP = .3 X 10 ^ W/Hz^ , a 100 nsec pulse, 

P 

a. 5 km visibility, a 1000 cm area receiver and K = .1 . For P = 0.5 and 
S/N = 1 for a single pulse we find the required energy vs range shown in 
Fig. 6 . Here a implies tuning exactly on resonance and implies 

tuning off resonance to satisfy Eq. (43), Here we have assumed that the natural 
CO abundance between 0.15 - 0,5 ppm (winter) and 0,3 to 3.0 ppm (summer )^^^ 




leads to an absorption of 0,96 Figure 6 clearly shows the advantage 

gained by optimizing the absorption cross section. Figure 6 assumes a CO 

concentration of 20 ppm over a length of 100 m in the path. For this case 

N.R < < N L up to a range of 10 km. Finally, Fig. 6 shows that for trans- 
A p 

mitted energies of 20 mj backscattered absorption measurements are possible 
for ranges up to 10 km. For most cases in urban areas, a 10 km range is more 
than adequate, and the single ended advantage of the topographical absorption 
measurement may offset the disadvantage of higher required laser power. 

An alternate way to compare the long path absorption measurement methods 
is to consider the minimum detectable pollutant concentration vs range. 

Figure 7 shows the result for the detection of SO^ , NO2 j and CO at the 
transmitted powers indicated. The methods have comparable sensitivities out 
to--a range where the S/N for the topographical detection begins to fa3J. due 
to lack of ^et^^rned signal. Here a S/N ratio of 100 is assumed which allows 
intensity changes to to be measured. The 100 mJ transmitted energy could 
be in a single pulse or equivalently in 100 averaged pulses of 10 mJ energy 
each. Figure 7 illustrates the very high sensitivity of the long path absorption 
methods and indicates the increase in transmitted power necessary to gain 
the advantage of a single ended absorption measurement. 

The long path absorption methods have been confirmed by experimental 
measurements,^^ Using diode lasers, Hinkley^'^^'^'^ 'has measured in 

automobile exhaust locally and over a parking lot area with a retroreflector 
return. Recently a diode laser system has been constructed for the monitoring 
of SOg- stack emissions. The operation of diode lasers at temperatures near 
77^K promises to extend considerably their applications to remote monitoring 
by absorption. 
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The topographical absorption method vas recently experimentally 
demonstrated by Henningsen et They used a 20 pj ^ 2.1 tim parametric 

oscillator source to remotely detect CO by topographical backscattering 
over a I 07 m range. Figure 8a shows the CO overtone spectrum taJcen in 
the laboratory with the LiNbO^ tunable parametric laser source. Figure 8b 
shows two rotation lines of CO taken over the 107 path. The measurement 
was not made differentially but by a continuous scan over a few minutes time. 
Kven so^ atmospheric scintillation did not reduce the. signal to noise signif- 
icantly, The measured CO detection sensitivity, s/lT , and range confirmed 
the theoretical results given by Eq. (4-2), Of course more extensive experi- 
mental studies need to be made, but the work to date has illustrated the 
useftillness of the long path absorption methods. 

E. Differential Absorption 

The differential absorption method overcomes the main disadvantage of 
the long path absorption methods - lack of depth resolution. By monitoring 
radiation backscattered from distributed Rayleigh and Mie particulate scattering, 
and tuning the transmitted wavelength on and off an absorption line of the 
pollutant molecule, depth resolved absorption measurements can be made with a 
depth resolution ct/2 as shown in Fig. 1. The differential absorption method 

thus combines in an optimum way the. depth resolution advantages of the Raman 
method with the very large absorption cross sections and resultant high 
sensitivity of the absorption methods. In this section the required transmitted 
energy, the signal to noise ratio and the sensitivity of the differential absorp- 
tion method is discussed. Finally, recent experimental results are reviewed to 
illustrate the advantages of this remote pollutant measurement method. 
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The differential absorption method was first suggested by Schotland 
lift 

in 1964'^ using a searchlight as a light source. The method was recently 
dlscuased by Igarashi^^ and first analyzed in detail by Eyer and Garbuny^^ 
and by Measures and Pilon-^ Wright^^ has also considered the application of 
differential absorption to pollution monitoring; In this section we follow 
the approach of Byer and Garbuny and first find the required transmitted 
energy, the signal to noise ratio and then the sensitivity of the measurement 
vs range. 

Before proceeding with the analysis, it is useful to schematically 

describe pollutant measurement by the differential absorption method. Figure 

9 illustrates the returned bachscattered power when the transmitted wavelength 

is tuned on and off the CO absorbing transition. The assumed transmitter 

end receiver characteristics are shown In the figure. The range derivative 

of equals the CO concentration shown at the bottom of the 

r ' r 

figure. The differential absorption method thus describes a differential in 
space with the differential tuning on and off the absorption line being under- 
stood. Figure 9 illustrates both the depth resolution and the sensitivity of 
the differential absorption method. 

To evaluate the required transmitted energy for a differential absorption 
measurement we proceed as previously for the topographical absorption method 
[cf Eqs. (37) - (^2)]. However, now the effective target reflectivity (p/rt) 
is replaced by ^ ct^(r) with total b.ackscatter coefficient given by 

3(R) = W 
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vhere and 3.^ are given by Eq. (it) and Eq. (6). However in the 

K Mlb 

infrared, caution must be used in applying Eq* 

* Table III shows that the effective reflectivity due to distributed 
particulate backscattering is approximately four orders of magnitude less 
than that for topographical targets. Thus the depth resolution advantage 
obtained by the differential absorption method is acquired at an increased 
required transsiitted energy. 

Substituting for p/jt in Eq. (42) we obtain for' the required-transmitted 
energy for the differential absorption measurement method . 




1 

op(R) 



exp 2[g^-^ f R 

1 - exp(-2 


(i^5) 


The’ required transmitted energy now varies as 1/vt instead ofvr as- for the 
topographical absorption method because of the bandwidth 'requirement for depth 
resolved measurements. Using Eq. (6) to find 3.,_(R) for a visibility of 
5 km and assuming the same parameters as previously for the topographical 
absorption of CO , we find the required transmitted energy shown in Fig. 10. 
Here is the CO cross section on resonance and ' the optimum 

cross section given by Eq. (43). At a given range, the differential absorption 
method requires about three orders of magnitude more transmitted energy than 
the topographical absorption method. A 10 mj pulse energy allows a depth 

g 

resolved absorption measurement to a 2 km range with a — = l,n = l and 
a l6 m depth resolution. The signal to noise ratio improves with increased 
transmitted energy but only as the square root of the number of averaged pulses. 
There is also a tradeoff in absorption sensitivity and depth resolution as 
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dl&cussed previously* Obviously a number of the variables can be adjusted 
to improve the range and sensitivity of the differential absorption method. 


However, Fig. 10 shows that transmitted energies in the 10 mJ - 100 mj range 
are necessary for differential absorption measurements in the infrared. The 
required energy is reduced somewhat in the visible due to increased scatter 
return and to more sensitive detectors. Examples of SO^ , NOg and Ig 
detection are given by Measures and Pxlon^ and the required energy for NOg 
detection is discussed by Byer and Garbuny.'^-^ 

The evaluation of the differential absorption method can be presented 
in an alternate way. If we assumed a fixed transmitted energy , we can 

ask how the signal to noise ratio varies as a function of range. For a dark 
ciirrent limited detector we "can invert Eq. (^5) solve for the S/lT ratio 
to find 



|(E) = Ej/E^(E) X « 1 
= 0AyEjR) 


(W) 


where x » 2 N L o and E^(R) is given , by Eq. (45) with S/N a 1. 

Figure 11 shows the signal to noise vs range for the detection of CO 
using an InSb detector and for the detection of KOg using a photomultiplier 
detector. Hereit is assumed that the maximum S/N = 100 which corresponds to 
a measurement accuracy 6 = Ax/x of 1^. For a 100 mJ transmitted pulse the 
Signal to noise ratio decreases beyond 1 km leading to a decreasing measurement 
sensitivity with increasing range. For example, for both CO and NOg the 
minimum measurable concentration for a 15 m depth resolution is 0.l4 ppm and 
0.85 ppm which holds to a 1 km range for a 100 mJ transmitted pulse. At 
2 km range the sensitivity is reduced to 1.4 ppm and 8.5 ppm. These sensitivities 
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can be improved by increasing the integration time and therefore the sampled 
depth and absorption length. However_, care must be taken to adequately 
suppress the background power. 

Schotland^^ has also considered errors involved in the differential . 

absorption measurement due to parameters other than detector signal to noise 

ratio. These parameters include atmospheric properties, uncertainties in the 

absorption coefficient and laser frequency, and signal to noise limitations of 

the returned power measurement. Wright has added digitization noise to the 

list of important contributions to the overall signal to noise ratio. However, 

122 

as pointed out by Schotland, it is ultimately the returned power signal to 
noise ratio which limits the range and accuracy of the differential absorption 
method. Within the past year experimental observations have verified the 
theoretically predicted parameters of the differential absorption method. 

118 

The earliest experimental work was performed by Schotland who used a 

temperature tuned Ruby laser to measure the water vapor vertical profile by 

the differential absorption method. Using a cw 3*39 HeNe laser Granatstein 

et al,,^^ demonstrated differential absorption measurement of CH|^ and COg 

in a controlled laboratory experiment.. The first reported differential 

119 

absorption measurements in the atmosphere were performed by Igarashi who 
used a dye laser to probe for NO^ over a 300m path. Concentrations of NOg 
and SOg were measured to 0.1 ppm over a 100m path using a 1 mj dye laser 
transmitter. 

Recently differential absorption measurements were considerably extended 

12k 125 

by Rothe, Brinkman and Walther ^ of Cologne, Germany. They used a 
1 mJ, 30 nsec grating tuned dye laser operating between k^'^Q - ^700 A 
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for the tunable laser transniitter. This spectral region was chosen to 
include the more pronounced NO^ absorption peaks. Following beam expansion 
the laser output was transmitted into the atmosphere and the backscattered signal 
was received with a 60 cm diameter Cassegrainian telescope and focused into 
a monochromator and onto a photomxaltiplier tube used in the photon counting 
mode. The differential range element was set at 5 OO m and ranges out to 
3.75 ^ were probed. Concentrations as low as 0.2 ppm were detected over 
the city of Cologne. 

The first calibration measurements of the differential absorption method 

12*5 

■were recently made by Grant et al, , Using a flashlamp pumped dye laser 

o 125 

operating near 4450 A at 4-8 mj pulse energy. Grant et al. , measured NO^ 

of known concentration in a sample chamber at 400 m range. Figure 12 shows 

the measurement results by the differential absorption method compared to the 

calibrated NOg concentration in the sample chamber. The uncertainty in the 

measurement is 20 ppm for a path length of 2.45 tn at S/N » 1. This is equivalent 

to 0.5 PPCQ in 100 m path. The error analysis for this measurement includes 

consideration of power signal to noise, atmospheric conditions and digital 

processing error. 

Perhaps the best demonstration of the differential absorption method 
is the two dimensional NO^ density plot over a factory in the city of Cologne 
Shown in Fig. 13* This composite of a series of range resolved measurements 

128 

made by Rothe- et al. , illustrate the potential of the differential absorption 
measurement for pollution monitoring over urban areas. For these measurements 
nearly 2000 pulses of 1 mJ each were averaged. The range resolution is 
about 100 m. For the plot shown a total of 40,000 pulses were used at 1 Hz 

repetition rate. 
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To date the differential absorption method has been experimentally 
verified only in the visible spectral range. However, plans are underway in 
a number of research laboratories to carry out measurements in the infrared. 

For example W^^lther, and Asai and Igarashi are planning to use a TEA- 
COg laser source for the detection of ozone near -10. [am. Also, work is 

129 

being completed on a high energy 1.4 pji to 4,4 pm parametric tunable source 
for application to near infrared differential absorption measurements. Based 
on previous experimental results in the visible, range resolved pollutant 
detection in the eye safe region of the infrared should prove to be a very 
useful monitoring tool in the near future. 

III. COKCLUSIOIT 

The advances in the technology of tunable laser sources has opened new 
possibilities for the remote measurement of atmospheric pollutants. From 
the analysis and examples presented in Section II it is now possible tO' determine 
the required transmitted energy for, each pollutant detection method and the 
sensitivity and range of the method at a given transmitted energy. 

The Raman method, due to the very small R«man cross section, requires 
a fixed frequency laser of approximately 1 J per pulse energy to achieve even 

a minimal range and sensitivity. The lack of sensitivity, required operation 
in the visible region where eye safety is a serious problem, and the background 
and interference signals preclude the use of the Raman method for trace pollutant 

detection. 

For local remote monitoring, the new four wave Raman mixing process which 
has scattering efficiencies 10 times that of spontaneous Raman scattering, 
looks very promising for trace gas analysis. However, experimental work 
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remains to be done to verify the sensitivity predictions. 

The resonance and fluorescence backscatter detection methods require 
a tunable laser source as do the absorption methods. From the analysis of 
these remote detection methods it is apparent that the resonance and fluor- 
escence backscatter method requires approximately 100 mj of transmitted 
tunable energy for detection in the infrared. This high required energy^ 
lack of depth resolution due to the long fluorescence decay times and variable 
cross section due to quenching, reduce the usefulness of the resonance method 
for remote monitoring in the troposphere. Only at high altitudes where quenching 
is not a problem and for atomic and molecular translsitons in the ultraviolet 
and visible does resonance backscatter offer detection advantages. For all 
other cases, the absorption methods are significantly better. 

The absorption methods considered are long path absorption using a retro- 
reflector or remote detector, long path absorption using a topographical target 
and the differential absorption method. For these methods the required power 
increases from approximately microwatts for the long path absorption method to 
pulsed energies of 1 mJ for the single ended topographical absorption method, 
to 10 to 100 mJ for the range resolved differential absorption method. For 
increased required tunable laser powers, the systems offer the advantages of 
single ended operation and single ended operation with depth resolution. The 
sensitivities of the three absorption methods are identical for equal path lengths. 

The widespx’ead use of the remote detection methods described in this 
paper depend critically upon the availability of narrow band, high energy 

130 

tunable laser sources in the spectral region of interest. Because of the 
molecular absorption spectra and the eye safety requirements the optimum 



spectral range is in the 2 pm to 15 ^ini infrared region. It is only recently 
that tunable sources with adequate pulse energy have been demonstrated in- 

131 

the infrared. However, with the development of high pressure gas lasers 

and high power nonlinear optical tunable sources, the prospect for remote 

air poUutioii measurements looks brighter. In addition, the demonstration 

l‘^3 

of diode laser operation in the infrared at liquid nitrogen temperatures 
significantly extends the possibilities for their widespread use for long 
path absorption meausrements and as heterodyne detector local oscillators 
to improve the sensitivity of the topographical and differential absorption 
methods. With the improvement of tunable infrared laser sources, remote air 
pollution measurement should find wider application in the near future. 
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figube captions 


1. Depth resolution for backscatter remote detection methods. The 
sampled depth is AR = ct/2 where t is the laser pulse width or 
detection integration time. 

2. Detectable CO concentration vs range by the Raman backscattering 
method. Here a doubled Ruby laser source is assumed. The other 
parameters are discussed in the text. 

o 

3. a) Raman return from ordinary air at a range of 30 m for a 337^ A 
transmitted beam from a pulsed nitrogen laser, 

b) Raman signal from an oil smoke showing the presence of , 

CI^ , HgS , CO , NO and SOg (after Kobaysi and Inaba, ref. 75). 

4. Fluorescence backscatter pollutant detection sensitivity vs range 
for CO in the infrared at 4.6 pm. A 100 mj transmitted pulse is 
assumed and a single pulse (n=l) at' S/N * 1 is evaluated. Sensitivities 
ere shown for optically thin and. optically thick pollutant layers, 

5. a) Calculated heterodyne detection sensitivities of SOg using a 
COg laser as a local oscillator vs the IF bandwidth- integration 
time product. 

b) Calculated heterodyne sensitivities for COg detection assuming lOC^ 
irradiation of the detector (solid line) and 20?5 irradiation of the detector 
(da shed line) and measured sensitivity (crosses) for 20 ^ detector 
irradiation. (After Menzies ref. 70). 
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6. Required transmitted energy vs range for CO detection at 4-. 7 

*^MAX i*^dicates timing on line center and indicates tuning to 

■ the optimum cross section for minimum required energy. A 20 ppa 
CO cloud 100 m thick is assumed in the path. (After Byer and 
Garbuny, reference 33). . 

7* Minimum measurable pollutant concentration vs range at s/N = 100 
for SOg f NOg and CO. The parameters are similar to those assumed 
in Fig. 6 except for a shot noise limited photomultiplier detector 
with F = Af = 20 MHz and 7] = for the SOg and NOg detection. 

8. a) Vibrational-rotational overtone, spectrum of CO by a tunable parametric 

laser. The measured cross section for CO at atmospheric pressure 
-PI P 

'V is 0.75 X 10 cm^. 

b) Two rotational absorption lines in CO of 18 atm-cm optical density 
at 107 m distance. Zero transmission is presented by the small horozontal 
lines at the beginning and end of the trace (after Henningsen et al. 
reference 115) • « ■> 

9. Backscattered power vs range for detection of CO by the differential 

absorption method. and P°^(R) are the returned powers when 

tuned off and on the CO absorption line. Also shown is the logarithm 
of the ratio (dashed line) and the range derivative of £n[P^^^/p°^] 
which equals the CO pollutant concentration. 

10. Required transmitted energy vs range f or . CO detection at 4.7 pm by 
the differential absorption method. For a 100 nsec resolvable pulse 
the depth resolution is I5 m and the minimum measurable pollutant 
concentration for S/N = 100 is 0.14 ppm. 
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11. S/N and measurement accuracy, 6, vs range for the differential 
absorption method assuming a 100 mJ, 100 nsec transmitted pulse. 

The depth resolution is 15 m. Here V is the visibility and t 

given by Eq. (^3), is the optimiim cross section. 

12. A calibrated differential absorption measurement of NO^ at 400 m. 

NOg measured by the differential absorption method vs NOg measured 

by transmissometer in the 2.4-5 ® sample chamber. These measurements 

o 

were taken at 44l8 and 4443 A (After Grant et al. , reference 125 ). 

13 . NOg distribution over a chemical factory as derived from differential 
absorption measurements at the indicated directions at an altitude of 

V 45 m . The concentratins are given in ppm (after Rothe et al. , 
reference 126 ). 
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